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Abstract  This  paper  presents  an  adaptive  thresholding  algorithm  that  can  be 
used  in  conjunction  with  the  multi-pass  GIP-based  editing  method  to  adaptively 
eliminate  ground  moving  target  returns  ("non-homogeneities")  from  the  training 
data  used  for  space-time  adaptive  processing  (STAP)  applications,  such  as 
adaptive  radars.  The  algorithm  exploits  a  property  of  the  generic  structure  of  the 
ordered  GIP  statistic,  the  origins  of  which  have  been  theoretically  developed  by  the 
authors,  and  a  single  adjustable  parameter  related  to  the  Type  I  error  of  incorrectly 
excising  target-free  training  data  to  adaptively  determine  the  thresholds  for  excising 
target-contaminated  training  data.  An  iterative  application  of  the  method  improves 
the  distinction  between  non-homogeneities  and  background  clutter,  which 
increases  the  number  of  target-contaminated  training  bins  excised.  Moreover,  the 
editing  method  has  been  extended  to  reduced  degree-of-freedom  (DoF)  STAP 
implementations  such  as  multi-bin  post-Doppler  STAP.  As  a  consequence,  the 
method  is  practically  implementable  for  realistic  scenarios  with  very  limited  a  priori 
knowledge  (i.e.,  only  the  number  of  DoFs  is  required). 

The  resulting  STAP  performance  improvement  with  the  new  technique  is 
demonstrated  by  applying  it  to  both  high-fidelity,  site-specific,  high-traffic  density 
simulation  data  with  a  realistic  radar  waveform  and  experimental  (MCARM)  data. 
Performance  metrics  include  SINR  Loss,  Pd  vs.  Pfa  curves,  and  association  of 
target  detections  with  road  locations  (for  the  experimental  data).  Results  with  the 
simulated  data  using  the  developed  GIP-editing  technique  essentially  recover  the 
target  detection  performance  of  clutter-only  training  data  when  an  ideal  waveform 
(i.e.,  no  range  sidelobes)  is  used.  The  presence  of  a  realistic  waveform,  as 
might  be  used  with  pulse  compression,  produces  range  sidelobes  that  degrade 
the  performance  of  the  editing  algorithm  as  a  function  of  range  sidelobe  level 
due  to  reduced  sensitivity.  However,  significant  improvement  in  target  detection 
performance  is  still  obtained  relative  to  no-editing.  Analysis  of  the  experimental 
MCARM  data  with  the  developed  GIP-editing  technique  indicates  that  the  number 
of  detections  is  about  doubled,  relative  to  the  unedited  results,  and  that  they  are  also 
better  associated  with  known  road  locations  in  the  scenario.  A  preliminary  analysis 
of  the  computational  complexity  of  the  editing  algorithm  is  also  performed. 
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Abstract  -  An  adaptive  thresholding  algorithm  is  presented 
that  can  be  used  in  conjunction  with  the  multi-pass  generalized 
inner  product  (GTP)-based  editing  method  to  eliminate  non¬ 
homogeneities  from  the  training  data  used  for  space-time  adap¬ 
tive  processing  (STAP)  applications  such  as  adaptive  radars.  The 
algorithm  exploits  a  property  of  the  generic  structure  of  the 
ordered  GIP  statistic  along  with  a  single  user-specified  parame¬ 
ter  related  to  the  Type  I  error  of  incorrectly  excising  target-free 
training  data  to  adaptively  determine  the  thresholds  for  excising 
target-contaminated  training  data.  The  performance  of  the 
method  is  demonstrated  using  both  high-fidelity  site-specific  sim¬ 
ulated  data  with  both  ideal  and  realistic  waveforms  as  well  as 
measured  data  from  the  Multi-Channel  Airborne  Radar  Mea¬ 
surement  (MCARM)  experiment. 

I.  Introduction 

The  deleterious  impact  of  non-homo geneous  training  data 
on  covariance  estimation  and  space -time  adaptive  processing 
(STAP)  [1]  performance  for  real-world  airborne  ground-mov¬ 
ing  target  indicator  (GMTI)  radar  data  has  been  widely 
reported  [2-6].  The  main  consequences  are  under-nulled  clut¬ 
ter,  unintentional  cancellation  of  targets  of  interest,  and 
reduced  detection  performance.  This  has  led  to  the  develop¬ 
ment  of  improved  training  data  selection  techniques  that  seek 
to  edit  out  bins  that  are  sufficiently  non-homogeneous  from 
the  bulk  of  the  training  data.  Such  non-homogeneity  detection 
(NHD)  methods  may  be  either  data  dependent  or  data  inde¬ 
pendent.  In  this  paper  we  focus  on  a  data-dependent  technique 
based  on  the  generalized  inner-product  (GIP)  [3-4]. 

The  GIP,  y ,  is  a  scalar  statistic  that  incorporates  both  power 
and  phase  information  of  the  data  [4].  It  is  defined  as 
y  =  x  R  x  where  x  is  the  MN  X  1  radar  data  (M  =  num¬ 
ber  of  pulses,  N  =  number  of  channels),  and  R  is  the  interfer¬ 
ence -plus-noise  covariance  matrix,  which,  in  practice,  must 
be  estimated  from  training  data.  “Large”  values  of  this  statis¬ 
tic,  equivalent  to  the  power  of  the  whitened  data  vector,  corre¬ 
spond  to  nonhomogeneities  (e.g.  ground-movers)  and  should 
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be  excised.  However,  the  presence  of  targets  in  the  training 
data  corrupts  the  covariance  estimate  in  the  GIP  as  it  does  the 
STAP  weight  calculation,  resulting  in  underprediction  of  the 
GIP  relative  to  its  true  value  (i.e.  using  the  ideal  covariance). 
This  increases  the  difficulty  of  distinguishing  between  targets 
and  clutter  [6].  To  address  this,  an  iterative  technique  (called 
the  “multi-pass”  algorithm  [3])  has  been  developed  whereby 
the  reduced  training  set  is  used  to  recompute  the  GIP  and  the 
process  of  excision  is  repeated  until  convergence  of  the  GIP 
statistic  as  a  function  of  range  is  achieved. 

Unfortunately,  it  is  not  clear  how  to  define  what  values  of 
GIP  are  too  “large,”  i.e.  what  are  the  thresholds  for  excision? 
A  survey  of  previous  work  reveals  two  basic  approaches. 
First,  a  prescribed  number  of  bins  may  be  removed  corre¬ 
sponding  to  the  outlier  values  of  GIP  [3,7].  Improvement 
results  but  it  is  not  clear  how  to  optimally  choose  the  number 
of  bins  to  remove.  A  second  approach  [7,  8]  is  to  use  statisti¬ 
cal  properties  of  the  GIP  derived  from  its  probability  distribu¬ 
tion  function  (pdf)  to  set  the  thresholds.  This  method  works 
well  if  the  assumed  pdf  accurately  represents  the  data.  How¬ 
ever,  the  presence  of  inhomogeneous  training  data  and  the  use 
of  techniques  to  compensate  for  the  absence  of  adequate  train¬ 
ing  data,  such  as  diagonal  loading,  can  significantly  alter  this 
pdf,  making  the  use  of  deterministic  limits  from  theoretical 
considerations  inaccurate  for  practical  use.  The  use  of  realistic 
waveforms,  which  produce  radar  return  range  sidelobes,  will 
tend  to  also  adversely  affect  this  editing  process. 

This  paper  presents  an  adaptive  thresholding  algorithm  that 
can  be  used  in  conjunction  with  the  multi-pass  GIP-based 
editing  method  [3]  to  eliminate  ground-moving  target  returns 
from  the  training  data  used  for  STAP  applications.  The  algo¬ 
rithm  does  not  require  knowledge  of  the  pdf  of  the  GIP,  nor  an 
a  priori  estimate  of  the  number  of  targets  to  be  removed.  The 
algorithm  exploits  a  property  of  the  generic  structure  of  the 
ordered  GIP  statistic  along  with  a  single  user-specified  param¬ 
eter  related  to  the  Type  I  error  of  incorrectly  excising  target- 
free  training  data  to  adaptively  determine  the  thresholds  for 
excising  target-contaminated  training  data.  The  performance 
of  the  method  is  demonstrated  using  both  high-fidelity  site- 
specific  simulated  data  as  well  as  measured  data  from  the 
MCARM  experiment  [9]. 

In  Section  II,  the  method  is  described  and  the  theoretical 
origins  of  the  GIP  properties  exploited  by  the  algorithm  are 


Sorted  GIP  Profile:  Ideal  Cov. 


presented.  Extensions  to  reduced  degree-of-freedom  (DoF) 
STAP  implementations  such  as  multi-bin  post-Doppler  STAP 
[10]  are  also  described.  Section  III  presents  results  using  site- 
specific  simulated  data.  Performance  metrics  include  SINR 
loss  and  P^  vs.  Pfa  curves.  In  Section  IV,  results  with  experi¬ 
mental  MCARM  data  are  reported.  Section  V  summarizes  the 
findings  and  outlines  areas  for  further  research. 

II.  Thresholding  Algorithm  and  Theory 

The  adaptive  thresholding  algorithm  is  illustrated  by  apply¬ 
ing  it  to  the  first  data  cube  released  under  DARPA’s  Knowl¬ 
edge-Aided  Sensor  Signal  Processing  and  Expert  Reasoning 
(KASSPER)  program  [11].  This  data,  from  the  KASSPER 
2002  Workshop  [12],  simulates  an  L-band  radar  with  parame¬ 
ters  similar  to  the  system  used  under  the  MCARM  program 
[9]  and  includes  site -specific  clutter  computed  using  digital 
terrain  elevation  data  (DTED).  Thus  this  data  set  represents  a 
generally  heterogeneous  clutter  environment.  Additional  real- 
world  elements  of  the  simulation  include  array  calibration 
errors  on  the  order  of  5-10  degrees,  clutter  discretes,  platform 
crab,  internal  clutter  motion  (ICM),  and  a  significant  number 
of  ground  moving  targets  in  the  radar  mainbeam.  The  simu¬ 
lated  system  uses  32  pulses  and  11  channels  for  352  total 
DoFs.  Details  of  the  simulation  technique  used  and  additional 
radar  parameters  are  provided  in  [12]. 

Figure  1  demonstrates  the  adaptive  thresholding  process 
with  the  GIP  statistic  computed  using  both  the  ideal  covari¬ 
ance  and  the  sample  covariance  (using  200  range  bins,  100  on 
either  side  of  the  bin  under  test,  and  0  dB  of  diagonal  loading 
relative  to  the  white  thermal  noise),  as  a  function  of  one-way 
radar  range  for  the  KASSPER  data.  Other  researchers  have 
previously  analyzed  the  rank-ordering  of  the  resulting  GIP 
statistic  [3,4].  We  observe  that  the  rank-ordered  GIP  has  a 
broad  “central”  region  that  is  nearly  linear  and  find  the  best- fit 
line  to  this  region.  The  upper  and  lower  thresholds  are  then 
determined  to  be  the  points  that  deviate  by  a  prescribed  rela¬ 
tive  fraction,  a  ,  from  this  straight  line.  Figure  1  illustrates 
that  while  the  GIP  profile  computed  using  the  sample  covari¬ 
ance  is  significantly  reduced  in  comparison  with  the  ideal 
covariance  profile  due  to  the  contaminating  effect  of  targets  in 
the  training  data,  the  sorted  profiles  have  very  similar  struc¬ 
ture.  Repeated  application  of  the  iterative  editing  process 
improves  the  number  of  target-contaminated  bins  that  are 
excised.  In  practice,  we  have  found  that  typically  at  most  five 
GIP  iterations  are  required  to  recover  the  performance  of  the 
sampled  result  using  clutter-only  data  and  an  ideal  waveform. 

A.  Ordered  Statistic  Theory 

Two  key  theoretical  issues  related  to  the  algorithm  are:  1) 
the  robustness  of  the  extended  linear  region  in  the  sorted  GIP 
graph,  and  2)  the  relationship  between  the  value  of  a  and  the 
properties  of  the  GIP  statistic. 

To  investigate  the  first  of  these  issues,  the  normalized  GIP 
statistic  was  examined,  y'  =  y/Q  =  (x  R  x)/Q ,  where 
Q  =  NM  is  the  total  number  of  DoFs.  In  the  following  we 
present  highlights  of  the  derivation  while  details  may  be 
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GIP  Range  Profile:  Sampled  Cov.  Sorted  GIP  Profile:  Sampled  Cov. 


Fig.  1.  Comparison  of  ideal  and  sampled  covariance  results 
for  the  adaptive  thresholding  algorithm. 

found  in  [15].  Assuming  target- free  Gaussian  data,  the  nor¬ 
malized  GIP  has  a  modified  gamma  pdf  [8]  and  a  resulting 
cumulative  distribution  function  (cdf)  F(  y )  that  may  be  com¬ 
puted  explicitly  but  is  not  expressible  in  closed  form. 

If  p  samples  of  the  normalized  GIP  are  acquired  and  sorted 
in  ascending  order,  the  pdf  of  the  rth  sample  (  1  <r<p  )  in 
this  ordered  list  (the  “rth  order  statistic”)  may  be  computed 
using  standard  techniques  [13], 

The  moments  of  the  rth  order  statistic,  particularly  the  mean 
and  variance,  are  of  interest.  The  qth  moment  can  be 
expressed  in  terms  of  the  inverse  of  the  cdf  F  (t)  called  the 
quantile  function 

oo  1 
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where  the  transformation  t  =  F(  y')  has  been  made  and  /p(0 
is  the  standard  beta  pdf  [14].  Analytical  evaluation  of  these 
moments  are  difficult  since  a  closed  form  expression  for  the 
quantile  function  of  the  normalized  gamma  pdf  does  not  exist. 
Instead,  a  Taylor  expansion  of  the  normalized  gamma  cdf  was 
formed  around  the  most  probable  value  of  y' , 
Y max  =  1  _  1/G  >  and  d  was  observed  that  the  second  order 
term  in  the  expansion  vanishes  if  Q  >  1  .  Thus,  around  the 
most  probable  value  of  y' ,  a  good  approximation  to  the  quan¬ 
tile  function  may  be  found  by  inverting  the  expansion  of  the 
cdf  to  first  order 
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where  P^iY  max)  ~  *lQ/ (2tt)  from  the  normalized  GIP  pdf 
and  F0(imax)  is  the  zeroth-order  term  in  the  Taylor  expan¬ 
sion  of  the  cdf  evaluated  at  the  most  probable  value  of  y' . 


Substituting  (2)  into  (1)  and  using  general  results  for  poly¬ 
nomial  moments  of  the  beta  pdf  [14]  allows  approximate  eval- 


uation  of  the  moments  of  the  rth  order  statistic.  For  example, 
an  approximation  to  the  mean  can  be  shown  to  be 


This  result  demonstrates  that  the  mean  of  the  rth  order  statistic 
of  the  normalized  GIP  around  the  most  probable  value  is  lin¬ 
early  related  to  the  order  statistic  number  r.  The  absence  of  the 
second  order  term  in  the  expansion  of  the  cdf  F(  y7)  explains 
the  broad,  central  linear  region  that  we  exploit  in  the  GIP  edit¬ 
ing  algorithm.  The  excellent  accuracy  and  robustness  of  these 
approximations,  for  both  the  linear  region  using  ]Lir  lin  and  for 
the  entire  domain  using  jll  r~F  (r/(p  +  1 ))  ,  has  been  dem¬ 
onstrated  with  the  KASSPER  data  using  the  ideal  covariance 
and  a  single  trial  of  data,  see  Figure  2.  We  note  that  this  simu¬ 
lated  data  is  non-Gaussian  due  to  the  implementation  of  the 
ICM  model  [12]  and  heterogeneous  terrain. 


Sorted  Normalized  GIP  Profile 


Fig.  2.  Comparison  of  sorted  GIP  computed  from  the 
KASSPER  data  to  theory  (3). 


The  thresholds  for  editing  are  determined  as  the  range  bins 
corresponding  to  the  order  statistics  that  deviate  from  the  lin¬ 
ear  approximation  by  a  given  fractional  error 
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It  is  assumed  that  any  ordered  statistic  that  exceeds  this  frac¬ 
tional  deviation  from  the  central  linear  region  (assumed  to 
contain  “homogeneous”  data)  contains  a  target  and  should  be 
excised.  In  [15],  guidance  for  choosing  a  is  provided  in 
terms  of  a  function  with  the  number  of  DoFs  and  desired  level 
of  statistical  Type  I  error  as  inputs.  Typical  values  of  a  range 
from  0.15  for  large  full-DoF  scenarios  (as  in  Figure  1)  to  0.25 
for  small  reduced-DoF  situations. 


B.  Reduced-DoF  STAP  Implementation 

For  many  systems  it  may  be  impractical  to  use  full- DoF 
STAP  for  reasons  of  computational  complexity  and  the 
absence  of  sufficient  sample  support  for  covariance  estima¬ 
tion.  In  such  situations,  it  is  desirable  to  reduce  the  number  of 
DoFs  to  a  manageable  number  while  minimizing  any  loss  in 
performance  relative  to  the  full-DoF  formulation.  The  adap¬ 
tive  thresholding  algorithm  may  be  straightforwardly  applied 
to  such  reduced-DoF  formulations  by  simply  computing  the 
GIP  with  the  reduced-DoF  variables  and  then  applying  the 
same  thresholding  procedure  to  the  resulting  GIP. 


An  example  of  a  reduced-DoF  technique  is  extended  fac¬ 
tored  or  multi-bin  post-Doppler  STAP  [10].  With  this  method, 
the  number  of  spatial  DoFs  are  preserved  while  the  temporal 
DoFs  are  reduced  to  a  specified  number  of  (typically  orthogo¬ 
nal)  Doppler  bins  for  each  cell  under  test.  The  thresholding 
algorithm  may  be  applied  to  the  GIP  results  for  each  Doppler 
bin  as  a  function  of  range  as  described  before.  As  a  conse¬ 
quence,  the  number  of  contaminating  targets  in  the  training 
data  is  now  a  function  of  Doppler  bin.  To  further  reduce  the 
potential  of  target  contamination,  a  “master  list”  for  each  bin 
is  compiled  as  the  intersection  of  allowed  range  bins  amongst 
the  multiple  Doppler  bins  used  for  the  DoF  reduction  for  the 
specified  cell  under  test. 

III.  Simulation  Results 

In  this  section,  the  performance  of  the  multi-pass  adaptive 
thresholding  algorithm  when  applied  to  the  previously 
described  high-fidelity  site-specific  KASSPER  simulated  L- 
band  datacube  is  analyzed  in  conjunction  with  both  full-DoF 
and  reduced-DoF  multi-bin  post-Doppler  STAP  implementa¬ 
tions.  The  full-DoF  (352  DoFs)  implementation  used  200 
local  samples  ( K  =  200  )  and  0  dB  of  diagonal  loading  rela¬ 
tive  to  the  white  thermal  noise  to  estimate  the  covariance.  Two 
reduced-DoF  implementations  were  examined  using  3  and  5 
post-Doppler  bins  (33  and  55  DoFs,  respectively).  A  training 
sample  size  of  K  =  2Q  range  bins  were  used  to  form  the 
sampled  covariance  in  both  cases  (i.e.  66  and  110  samples, 
respectively)  with  0  dB  of  diagonal  loading.  The  same  train¬ 
ing  strategies  were  used  for  both  the  GIP  pre-processing  and 
STAP  steps.  The  value  of  the  thresholding  fraction  was  cho¬ 
sen  to  be  a  =  0.15  for  the  full-DoF  case  and  increased  to 
a  =  0.20  for  the  reduced-DoF  cases.  The  GIP  processing 
step  was  repeated  until  the  training  list  had  converged  for  each 
case  (approximately  5  iterations). 

Elsewhere  it  has  been  shown  that  the  GIP-editing  process  is 
effective  in  removing  high  SNR  targets  from  the  training  data 
for  both  full-  and  reduced-DoF  STAP  implementations  with 
this  data  [15].  A  consequence  of  this  editing  step  is  an 
improvement  in  the  ratio  of  signal-to -interference-plus-noise 
to  signal-to-noise  (SINR  loss).  This  is  demonstrated  in  Figure 
3  where  the  SINR  loss  as  a  function  of  range  and  Doppler  for 
all  three  cases  are  shown.  Note  that  the  mainbeam  clutter 
notch  is  located  at  the  two-way  Doppler  shift  of  ~51  m/s  due 
to  the  antenna  beam  being  steered  away  from  broadside.  The 
improvement  in  SINR  loss  with  editing  is  significant,  particu¬ 
larly  around  the  highly  contaminated  region  of  0  m/s  Doppler, 
where  the  improvement  with  editing  may  be  as  large  as 
22.5  dB.  Other  Dopplers  where  improvement  is  seen  include 
90  m/s  and  30  m/s  (near  to  the  clutter). 

To  further  investigate  the  benefit  of  the  GIP-editing  pro¬ 
cess,  the  KASSPER  space-time  beamformed  clutter  data  was 
post  processed  with  a  median  constant  false  alarm  (CFAR) 
algorithm  [16]  that  employed  100  range  bins  and  a  single 
Doppler  bin  to  compute  a  background  noise  level  for  the 
beamformer  output.  The  normalized  beamformer  output  was 
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Fig.  3.  SINR  loss  for  unedited  (top)  and  GIP-edited  (bottom) 
training  data.  Left:  full-DoF;  Middle:  5 -Bin  post-Doppler; 
Right:  3 -Bin  post-Doppler.  Two-way  Doppler  is  shown. 

thresholded  and  the  number  of  false  alarms  was  recorded.  The 
sample  processing  steps  were  then  applied  to  a  series  of  test 
targets  in  all  range  bins  for  a  specified  Doppler  shift.  The  tar¬ 
get  SNR  was  set  to  25  dB  at  the  closest  range  bin  which  repre¬ 
sents  a  radar  cross  section  of  approximately  5  dBsm.  The 
target  SNR  was  then  reduced  versus  range  using  an  inverse 
range  to  the  fourth  power  rule  which  resulted  in  a  target  SNR 
of  1 9  dB  at  the  longest  range  in  the  simulation.  The  number  of 
targets  exceeding  the  thresholds  when  processed  using  the 
same  beamformer  weights  and  CFAR  normalization  was  also 
recorded.  vs.  Pfa  curves  were  then  generated  using  the 
observed  number  of  false  alarms  and  detections. 

Figure  4  shows  the  P^  vs.  Pfa  curves  for  a  Doppler  velocity 
of  0  m/s  (note:  clutter  at  51  m/s)  around  which  many  targets 
may  be  found  (see  Figure  3).  The  unedited  results  for  all  three 
sampled  cases  show  dramatic  degradation  relative  to  the  opti¬ 
mal  ideal  covariance  case  due  to  the  impact  of  targets  in  the 
training  data.  However,  the  GIP-edited  results  demonstrate 
significant  improvement  of  Pd  relative  to  the  unedited  results 
due  to  the  effective  removal  of  contaminating  targets  from  the 
training  data. 

The  results  presented  so  far  have  used  an  ideal  waveform 
(i.e.  a  rectangular  pulse  whose  length  equals  the  inverse  of  the 
system  bandwidth  and  has  infinitely  low  sidelobes).  However, 
a  realistic  waveform  results  in  energy  leakage  (i.e.  sidelobes) 
of  the  target  response  in  range.  An  example  of  this  is  the  com¬ 
pressed  linear  frequency  modulated  (LFM)  waveform,  which 
was  used  in  this  study  and  is  described  below.  The  resulting 
increased  self-contamination  of  the  covariance  estimate 
reduces  the  value  of  the  GIP  making  it  even  more  difficult  to 
distinguish  non-homogeneities  from  the  background.  This 
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Fig.  4.  P^  vs.  Pfa  curves  for  Doppler  velocity  0  m/s.  Unedit¬ 
ed  (dashed)  and  GIP-edited  (solid)  training  data  results  for  all 
three  cases  shown.  Full  DoF  results  (FD)  use  K  =  200 , 5- 
bin  (PD5)  results  use  K  =  110  ,  and  3 -bin  (PD3)  K  =  66  . 

desensitization  of  the  GIP  causes  fewer  target  contaminated 
bins  to  be  edited  out  and  results  in  degraded  performance  rela¬ 
tive  to  the  edited  ideal  waveform  case. 

This  effect  was  investigated  using  the  KASSPER  data  with 
the  targets  reproduced  using  a  compressed  LFM  waveform 
with  a  Taylor  pulse  compression  taper  (with  nbar=6)  and  a 
duty  factor  of  0. 1.  Taper  sidelobe  levels  of  20,  25,  30,  35,  and 
40  dB  below  the  peak  response  were  examined.  The  previ¬ 
ously  described  5-bin  post-Doppler  (K=110)  STAP  architec¬ 
ture  with  four  guard  bins  on  either  side  of  the  cell  under  test 
was  used.  As  expected  the  GIP  variation  and  level  with  the 
converged  list  was  seen  to  increase  with  reduced  sidelobe 
level.  Pd  vs.  Pfa  performance  results  using  the  median  CFAR 
and  the  previously  described  target  inject  process  for  all 
waveforms  including  ideal  are  shown  in  Figure  5.  In  these 
results,  the  clutter  contribution  is  processed  with  the  ideal 
waveform  for  all  cases  while  the  LFM  waveform  affects  only 
the  targets.  Both  the  unedited  and  GIP-edited  results,  pre¬ 
sented  separately,  are  seen  to  degrade  with  increased  sidelobe 
level.  However,  except  at  the  largest  sidelobe  (20  dB)  and 
lowest  Pfa  levels,  GIP-editing  provides  a  significant  improve¬ 
ment  in  Pd  over  no-editing  with  the  LFM  waveform.  We  also 
note  the  insensitivity  to  further  improvement  as  the  sidelobe 
level  is  reduced  below  35  dB. 

A  potential  reason  for  the  loss  of  performance  with  GIP- 
editing  at  low  levels  of  Pfa  is  that  the  same  level  of  sample 
support  was  used  for  all  cases  resulting  in  some  loss  of  train¬ 
ing  data  locality  for  the  edited  cases  and  a  resulting  increased 
sampling  error  for  the  covariance  estimate.  This  was  con¬ 
firmed  by  performing  the  same  analysis  using  uncontaminated 
(i.e.  clutter-plus-noise-only)  training  data  without  and  with 
the  GIP-edited  training  list  and  the  same  training  strategy.  The 
resulting  Pfa  curves  for  both  were  very  similar  to  those  using 
contaminated  training  data.  This  indicates  that  the  loss  in  Pfa 
performance  with  the  edited  training  list  is  primarily  due  to 
non-locality  effects  in  the  training  data  rather  than  residual 
target  contamination  effects.  In  the  future,  we  will  examine 
the  degree  to  which  this  effect  is  site-specific  (i.e.  pronounced 
in  areas  of  significant  terrain  heterogeneity)  and  whether  it 
may  be  mitigated  by  using  sample  support  reduction  tech¬ 
niques  (e.g.  [22]). 
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Fig.  5.  P^  vs.  Pfa  curves  for  Doppler  velocity  0  m/s  with  PD 5. 
Unedited  and  GIP-edited  results  for  all  LFM  waveforms 
shown.  Legends  indicate  waveform  type/pealc  side  lobe  level. 


IV.  MCARM  Experimental  Results 

We  now  briefly  summarize  results  of  applying  the  GIP -edit¬ 
ing  method  to  the  data  from  MCARM  flight  5,  Acquisition 
575  [17].  Radar  parameters  for  this  L-band  scenario  may  be 
found  in  [9].  The  reduced-DoF  3-bin  post-Doppler  STAP 
implementation  in  conjunction  with  all  22  spatial  elements 
was  used  to  process  the  data.  The  covariance  was  estimated 
using  twice  the  number  of  DoFs  (i.e.  132  local  training  bins), 
0  dB  diagonal  loading  consistent  with  the  estimated  thermal 
noise  level,  and  2  guard  bins.  The  data  was  pulse  compressed 
using  a  1  MHz  LFM  waveform. 

Examination  of  the  resulting  rank-ordered  GIP  profiles  at 
various  Doppler  velocities  revealed  the  same  central  linear 
structure  as  found  in  the  simulated  data  and  predicted  from 
theory.  The  multi-pass  editing  algorithm  was  applied  with  the 
fractional  threshold  parameter  set  to  a  =  0.25  . 

The  resulting  GIP-edited  training  list,  along  with  the 
unedited  list,  was  used  to  perform  STAP.  Elsewhere,  the 
improvement  in  SINR  loss  profiles  as  a  function  of  Doppler 
velocity  with  GIP-editing  over  unedited  training  data  was 
illustrated  [15].  Median  CFAR  processed  beamformed  outputs 
at  two  Dopplers  (-40.3  and  -32.3  m/s;  clutter  at  -12  m/s) 
where  significant  improvement  in  SINR  loss  was  previously 
found  are  presented  in  Figure  6.  A  threshold  of  15  dB  was 
used  to  denote  detections.  At  both  Dopplers,  the  GIP-edited 
results  indicate  detections  that  were  not  present  in  the 
unedited  results.  Although  ground  truth  for  the  targets  in  the 
scenario  was  not  available  for  comparison,  we  do  note  that 
these  detections  correspond  in  range  to  the  indicated  roads. 

This  median  CFAR  detection  process,  with  a  threshold  of 
15  dB,  was  applied  to  all  processed  Doppler  bins  and  potential 
target  detections  were  identified.  Due  to  the  absence  of 
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Fig.  6.  Median  CFAR  processed  beamformed  output  for  the 
MCARM  experimental  data  with  both  unedited  and  GIP-ed¬ 
ited  training  data.  A  detection  threshold  of  15  dB  and  names 
of  known  roads  close  to  the  indicated  detections  also  shown. 

ground  truth,  the  detections  were  further  processed  from 
range-Doppler  to  range-azimuth  and  ultimately  latitude -longi¬ 
tude  using  the  WGS84  Earth  model  [18].  They  were  then 
compared  with  locations  of  known  roads  (using  TIGER/Line 
[19]  data).  It  is  reasonable  to  assume  that  detections  located 
by  this  process  that  are  close  to  roads  are  likely  actual  ground 
movers.  Figure  7  presents  the  results  of  this  analysis  for  both 
the  unedited  and  GIP-edited  training  lists.  First  of  all,  the 
number  of  detections  effectively  doubled  (from  69  to  134) 
with  the  use  of  GIP-editing.  Next  we  note  that  the  GIP-edited 
results  shows  greater  concentration  and  density  of  targets 
around  the  roads  than  does  the  unedited  results.  We  conclude 
that  GIP-editing  has  improved  performance  for  this  scenario. 

V.  Conclusions 

This  paper  presents  an  adaptive  thresholding  algorithm  that 
can  be  straightforwardly  used  in  conjunction  with  the  multi¬ 
pass  GIP-based  editing  method  to  eliminate  non-homogene¬ 
ities  from  the  training  data  used  for  STAP  applications.  Imple¬ 
mentations  for  both  full-  and  reduced-DoF  processing  were 
presented.  The  algorithm  exploits  the  fact  that  the  rank- 
ordered  GIP  statistic  has  a  broad,  central  linear  region,  the 
theoretical  origin  of  which  was  described  in  this  paper.  The 
algorithm  requires  a  single  input  parameter  a  the  allowable 
fractional  deviation  from  this  linear  region,  to  determine  the 
thresholds  for  excision.  Guidance  for  choosing  this  parameter 
is  provided  in  terms  of  a  function  with  the  number  of  DoFs 
and  desired  level  of  statistical  Type  I  error  used  as  inputs. 

Robustness  and  performance  of  the  method  was  demon¬ 
strated  by  applying  it  to  STAP  implementations  with  varying 
number  of  DoFs  and  sample  support  as  applied  to  both  high- 
fidelity  site-specific  simulated  data  and  experimental  data. 


MCARM  Data:  SMI  GIP  Edit 


Fig.  7.  MCARM  detection  results  (x’s)  in  relation  to  roads  in 
the  scenario.  Top:  No  edit.  Bottom:  GIP  edit.  Detections  in¬ 
creased  from  69  to  134  with  editing.  Mainbeam  of  radar  also 
shown. 

STAP  with  GIP-edited  training  data  was  shown  to  have 
improved  detection  performance  over  unedited  results  for  all 
cases.  Editing  performance  was  shown  to  degrade  in  the  pres¬ 
ence  of  LFM  waveforms  with  increasing  sidelobe  levels  but 
improvement  over  no -editing  was  still  significant. 

A  preliminary  analysis  of  computational  complexity  indi¬ 
cates  that  the  most  efficient  way  to  implement  the  iterative 
editing  process  when  using  a  local  “sliding  window”  strategy 

[20]  is  to  use  a  Choleslcy  factor  up/down-dating  scheme  in  the 
data  domain,  such  as  described  in  [21].  We  note  that  if  suffi¬ 
cient  storage  is  available,  the  resulting  converged  Choleslcy 
decomposition  for  each  range  bin  from  the  GIP-editing  pro¬ 
cess  may  be  used  in  the  subsequent  data  domain  STAP  step, 
thereby  improving  its  efficiency. 

Finally,  since  the  number  of  contaminating  targets  in  the 
training  data  naturally  is  reduced  when  lower  sample  support 
is  used,  we  are  also  investigating  ways  of  combining  the  GIP 
editing  approach  with  methods  that  minimize  sample  support 
requirements  such  as  knowledge- aided  pre -whitening  tech¬ 
niques  [22]. 
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Problem  Illustration 


*  Ground  moving  detection  performance  with  STAP  is  degraded  due  to 
inaccuracy  of  noise-plus-interference  covariance  estimate 


•  Two  common  sources  of  error  due  to  use  of  training  data  are  i)  sample 
error,  ii)  target  contamination 
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Sample  error  widens  clutter  notch 
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Significantly  degraded  Pfa  vs.  Pd 
performance  away  from  clutter 
may  result 


ISL 


ASAP2004/CMT  -03/  04  -  3 


Solution  Approach  Overview 


training  data 
selection 


Statistical 
_  Test  based  on 
Available 

Radar  data  (e.g.,  GIP) 


•  The  goal  is  to  remove  training  samples  that  contain  ground  targets 
from  the  training  set 

•  Use  multi-pass  Generalized  Inner-Product  (GIP)  test  with  adaptive 
thresholding  to  determine  outliers  which  are  then  excised  from  the 
training  data 

•  Use  resulting  edited  training  data  to  estimate  noise+interference 
covariance  matrix  as  required  for  adaptive  clutter  filter  application 

•  Improved  detection  performance  results 
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Objectives 

•  Develop  a  practically  implementable  adaptive  thresholding 
algorithm  for  use  with  the  multi-pass  GIP-based  editing 
method 


•  Establish  theoretical  understanding  and  robustness  of  key 
features  of  algorithm 


•  Extend  to  reduced-DoF  implementations 


•  Demonstrate  performance  improvement  using  GIP-based 
editing  algorithms  with  STAP  for  broad  range  of  operating 
environments  and  radar  systems 
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Outline 


Introduction 

GIP-Based  Editing 

-  Definition 

-  Illustration  of  Algorithm 

-  Theory 

Results 

Summary  &  Future  Work 
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Data-Dependent  GIP-Based  Editing 


The  Generalized  Inner  Product  (GIP)  statistic,  y,  may  be  used 
to  identify  non-homogeneities  in  the  data 


2 

x  =  data 

/(k)  =  xH  (k)R~\k)x(k)  = 

R“1/2x 

?  k  =  range  bin 

R  =  clutter  cov. 

Outlier  value  of  y  indicates  that  clutter-whitened  data 
differs  from  noise;  non-homogeneity  present 


Works  well  when  covariance  estimate  accurate  (i.e.  ideal 
covariance);  distinction  between  noise  and  targets  clear 


When  GM  contaminated  training  data  used  to  estimate  R, 
value  of  outlier  y  reduced,  makes  detection  much  more 
difficult  since  distinction  between  noise/targets  not  as  clear 
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How  to  Determine  Thresholds  for  Excision? 

•  Remove  prescribed  number  of  bins;  fix  thresholds  empirically 

-  How  to  generalize/optimize  ?? 

•  Deterministic  approaches  based  on  ideal  statistical  properties 

-  May  be  significantly  in  error  due  to  practicality  issues  (e.g.  non¬ 
stationary  training  data,  loading  of  covariance  estimate,  etc.) 


Our  Approach: 


•  Analyze  sorted  GIP  statistic,  edit,  iterate 

-  “Multi-pass  GIP-editing”;  [Melvin  etal.,  RadarCon  96/97]  and  others 

•  Apply  adaptive  thresholding  to  determine  and  excise 
outliers  from  training  data 

-  Exploits  characteristic  structure  of  sorted  GIP  statistic 

•  Iteration  improves  performance  by  enhancing  distinction; 
typically  converges  <  -  5  iterations 

•  Has  been  extended  to  reduced-DoF  (e.g.  multi-bin  post- 
Doppler)  STAP 

-  Perform  editing  using  transformed  data  on  per  Doppler  bin  basis 
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Order-Based  GIP  Algorithm 


Ideal  R  result 

Large  y  change 
when  GM 
present 

Linear  Fit  to 
central  region  of 
sorted  y 

Adaptive  limits 
shown  (due  to 
deviation  from 
linear  fit) 

a  =  0.15 


Sampled  w/GM 
result 

Small  y  change 
when  GM 
present 

Same  sorted 
structure  as 
with  Ideal  R 

Adaptive  limits 
shown  (with 
same  limit 
parameter,  a) 


GIP  Range  Profile:  Ideal  Cov. 


Sort 


range  (km) 

GIP  Range  Profile:  Sampled  Cov. 


Sort 


:  L-band  Example 


Sorted  GIP  Profile:  Ideal  Cov. 


Sorted  GIP  Profile:  Sampled  Cov. 
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Summary  of  Theory  of  Rank-Ordered  GIP 

Existence  of  broad  central  linear  region  a  consequence  of  the 
properties  of  the  order  statistics  of  the  GIP 

-  The  PDF  of  the  rth  ordered  sample,  out  of  p  total,  is  well  known** 


The  moments  of  the  rth  order  statistic  can  be  expressed  in 
terms  of  the  inverse  of  the  GIP  CDF,  F,  e.g.  for  the  mean 


Beta  PDF 


Assuming  Gaussian  data,  exact  expressions  not  available  but 
possible  to  show  that  second  order  term  in  Taylor  expansion 
of  CDF  around  most  probable  value  vanishes  so  that 


Mr~F 


-1 


\2k 

r  \ 

+  1- 


Q  =  # DoFs 


Thus  the  mean  of  the  rth  ordered  GIP  statistic  well 
approximated  as  linear  function  of  normalized  rank  order 
around  most  probable  value 


H.A.  David,  Order  Statistics,  1981. 
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norm.  GIP  (y/Q) 


Comparison  of  Theory  to  Simulated  Data 


•  Theory  compared  to  single  trial  of  GIP  computed  from  L-band 
KASSPER  data  using  the  ideal  covariance 

-  In  general  non-Gaussian  due  to  Internal  Clutter  Motion  (ICM) 
model,  heterogeneous  terrain 

Sorted  Normalized  GIP  Profile 


•  Extent  of  linear  region  and  agreement  with  theory  despite 
non-Gaussian  nature  of  data  indicates  theory  is  somewhat 
robust 
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Choice  of  Threshold  Parameter,  a 

GIP  analysis  suggests  choice  of  a  for  given  number  of  DoFs  to 
achieve  a  desired  Type  I  Error  probability  (i.e.  probability  of 
mistakenly  editing  out  a  range  bin  that  is  target-free) 


For  full-DoF  simulations  with  #  DoFs  >100  typically,  a  =  0.15  works  well 

For  reduced  DoF  (e.g.  post-Doppler)  should  increase  value  of  a  to  allow 
for  greater  variance  of  GIP  (“fatter”  tails) 

Use  as  guide  for  scenarios  with  non-stationary  data 

ISL 
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Outline 


Introduction 

GIP-Based  Editing 

Results 

-  High-Fidelity  Site-Specific  Simulated  Data  (KASSPER#1) 

-  Experimental  Data  (MCARM) 

Summary  &  Future  Work 
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First  DARPA/AFRL  KASSPER  Scenario 


Rx:  36.27°  N,  117.8°W 
3000  m,  (2088  m) 
Tx:  36.27°  N,  117.8°W 
3000  m,  (2088  m) 
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INFORMATION 
SYSTEMS 
LABORATORIES,  INC. 


www.islinc.com 


0  10  20  30  40  50  60  70  80  90  . 

(35.82°, 118.3°)  x  distance  (km) ->  (longitude)  2i-Mar-2002 


(36.72°  N,  117.2°  W) 


Mainbeam 
steering  direction 
shown 

BW=10  MHz, 
PRF=1984  Hz, 

11  subarray 
elements, 

32  pulses;  L-band 

Platform  is 
heading  west  at 
100  m/s  and  3km 
altitude 

Range  swath  of 
interest  shown 

About  200  targets 
in  the  mainbeam 

Array  calibration 
errors,  discretes 
included 

3  degree  crab 

15  mph  average 
wind  speed  (ICM) 
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Full-DoF  Results  :  SINR  Loss 


Doppler  (m/s) 

Q  =  352,  SMI  with  200 
training  bins,  0  dB  dl 

Single  GIP  editing  pass 
shows  only  marginal 
improvement 

Converged  sampled  result 
very  similar  to  sampled 
result  without  ground 
movers  in  the  covariance 
estimate 

GIP  results  also  indicate 
that  ground  movers  have 
been  very  effectively 
edited  out  of  training  data 


Doppler  (m/s) 


Doppler  (m/s) 


SMI200  w  GM,  Single  Pass  GIP  EditSMI200  w  GM,  Converged  GIP  Edit 


-100  0  100  -100  0  100 
Doppler  (m/s)  Doppler  (m/s) 


ASAP2004/CMT  -03/  04  -  1 5 


SINR/SNR  (dB)  SINR/SNR  (dB) 


no.  of  contaminating  targets  no.  of  contaminating  targets 


Target  Contaminated  Range  Bin 

Full-DoF 

No.  of  Targets  in  Training  Data 


10 


20  30  40  50 

target  SNR  (dB) 

3,5-Bin  Post-Doppler 

No.  of  Targets  in  Training  Data 


60 


10 


20  30  40  50 

target  SNR  (dB) 


60 


Removal 

Extent  of  removal  of 
target  contaminated 
range  bins  from 
training  data  with  GIP 
editing  shown 

Scenario  contains  206 
detectable  targets 
(SNR  >  12  dB)  within 
169  range  bins 

For  Full  DoF,  20-80 
contaminating  targets 
per  range  bin  pre¬ 
editing;  <18  after  GIP 
editing 

Remaining  targets  tend 
to  have  low  SNR 

Same  trend  with  post- 
Doppler  but  fewer 
contaminating  targets 
due  to  reduced  sample 
support 

Very  few  low  SNR 
targets  remain  after 
editing  III 
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CFAR  Detections:  SMI  without  Editing 


CFAR  Processed 
Beamformed  Output 


Detection  Information 


-100 


0  100 
Doppler  (m/s) 


17  Associated  Detections,  17  False  Alarms 


•  5-Bin  post- 
Doppler  STAP, 
K=110  +  0  dB  dl 


•  Distortionless 
weights  with 
median  CFAR 
target  detection 

•  Threshold  =  13 
dB;  single 
detection  per 
range  bin 

•  Detected  targets 
associated  with 
ground  truth 
when  possible 

-  40  m  in  range 

-  8  m/s  in 
Doppler 


Poor  target 
detection  due  to 
target 

contamination 
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CFAR  Detections:  SMI  with  GIP  Editing 


CFAR  Processed 
Beamformed  Output 


-100  0  100 


Doppler  (m/s) 


Detection  Information 


38  - 


36  - 


M 

3  ® 


X  Truth 

•  Raw  Detections 
O  Associated  Detections 
x  Associated  Truth 


X 

*0 

s 

_ jp 

© 

1 

r 

4  : 
© 

i 

-50  0  50  100 

Doppler  Velocity  (m/s) 


143  Associated  Detections,  22  False  Alarms 


•  Same  STAP 
processing, 
CFAR  detection 
process  as  in 
previous  slide 

•  Dramatic 
improvement  in 
number  of 
associated 
detections  due 
to  GIP  editing 

•  Comparable 
false  alarm  rate 
as  with  unedited 
results 

•  Clutter-Only  SMI 
performance: 
147  AD/ 13  FA 
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Pd  vs  Pfa  Performance  :  High  Velocity  Target 


Pf  vs.  Pd  :  SNRq  =  25  dB,  Doppler  =  0  m/s 


1 000  injected  test  targets,  all  ranges,  at  Doppler  velocity  =  0  m/s  no  array  cai 

(clutter  at  51  m/s),  target  SNR  is  25  dB  at  closest  range  bin  (~  5  dBsm) 

Detector  includes  median  CFAR  normalization  of  beamformed  output  prior 
to  thresholding 

“No  Targ”  =  no  targets  in  training  data  (clutter-only); 

“No  Edit”  =  no  editing  of  target  contaminated  training  data; 

“GIP  Edit”  =  GIP  edited  training  data 

No  targets  in  CFAR  training  data  for  all  cases 


GIP  Editing  recovers  performance  to  close  to  clutter-only  results 
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Pd  vs  Pfa  Performance  :  Low  Velocity  Target 

P  vs.  Pd  :  SNRo  =  25  dB,  Doppler  =  24.9021  m/s 

—  Ideal  Cov. 

""  FD,K=200,  No  Targ 

-  -  FD,K=200,No  Edit 

-  FD,K=200,GIP  Edit 

"  PD5,K=1 10,  No  Targ 

-  -  PD5,K=1 10,  No  Edit 
-  PD5,K=1 10,GIP  Edit 

PD3,K=66,  No  Targ 

-  -  PD3,K=66,  No  Edit 

-  PD3,K=66,GIP  Edit 

—  Model  Only 

-  -  Conv.  Doppler _ 


-4  _?  _i 

10  10  10  10 

Pfa 

•  Same  analysis  as  previous  slide  but  targets  now  injected  at  Doppler 
velocity  24.9  m/s,  which  is  closer  to  the  clutter  ridge  (51  m/s) 

•  “No  Targ”  results  show  degradation  due  to  sample  estimation  errors 
while  “No  Edit”  results  are  further  degraded  due  to  target  contamination 

•  “GIP  Edit”  results  approach  “No  Targ”  (clutter-only)  performance  as 
tolerable  Pfa  increases 

•  Loss  of  “GIP  Edit”  performance  at  low  Pfa  a  result  of  wider  spreading  of 
training  data  over  range  leading  to  higher  sample  estimation  errors  due 
to  terrain  heterogeneity 

•  Examine  if  can  be  improved  by  combining  with  techniques  that  reduce 
sample  support  requirements  (e.g.  knowledge-aided  pre-whitening) 

ASAP2004/CMT  -03/  04  -  20 


Experimental  MCARM  Scenario 


Rx:  39.38°  N,  75.97°  W 
3071.1888  m,  (1  m) 
Tx:  39.38°  N,  75.97°  W 
3071.1888  m,  (1  m)60 
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(40.14°  N,  75°  W) 


0  20  40  60  80  100  120  140  160 

(38. 62°, 76. 95°)  x  distance  (km)  ->  (longitude) 


O  5 

(38.62a,76.95°) 


1 5  20  25  30 

x  distance  (km)  — >  (longitude) 


Multi-Channel  Airborne 
Radar  Measurement 
(MCARM)  scenario  (Rome 
AFB  RL-TR-96-49) 

Analyze  MCARM  Flight  5, 
Acquisition  575 

L-Band  (1240  MHz) 

East  of  Baltimore,  platform 
looking  easterly 

PRF=1984  Hz,  2x11 
channels,  128  pulses 

Many  major  roads 
intersecting  mainbeam 

Low  Bandwidth=0.8  MHz 
(~190  m  range  resolution, 

346  bins) 

Typical  GMTI  radars  use 
higher  BW,  range  resolution 

Calibrated  steering  vectors 
provided 

Data  is  pulse  compressed 
using  1  MHz  LFM  waveform 

ISL 
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Experimental  MCARM  Data  SINR  Loss 


SINR  Loss,  Range  =  21 .88  km 


SINR  Loss,  Range  =  30.88  km 


no  edit 
GIP  edit 


-25  0  25 

two-way  Doppler  (m/s) 


-25  0  25 

two-way  Doppler  (m/s) 


50 


Hwy 

13 


•  SINR  Loss  using  3-Bin  post-Doppler  with  2xDoF  sample  support 
(132  training  bins)  with  0  dB  diagonal  loading  and  2  guard  bins 

•  Results  at  ranges  consistent  with  known  roads 

•  Unedited  and  GIP-edited  results  shown 


•  Improvement  with  GIP  approaches  13  dB  for  some  Dopplers  at 
shorter  range,  more  moderate  at  longer  range 

•  GIP-edit  results  indicate  detections  at  ranges  corresponding  to 
known  road  locations  not  present  with  unedited  results 
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y  distance  (km)  ->  Latitude  y  distance  (km)  ->  Latitude 


MCARM  Experimental  Detection  Results 


mcarm  Data:  smi  No  Edit  69  Detections 


20  25  30  35  40  45  50  55  60  65  70 


x  distance  (km)  ->  Longitude 


No  ground  truth 

-  Can’t  perform  target 
association 

Instead,  find  max 
response  in  azimuth 
for  each  target 

Convert  from  range- 
az  to  lat-lon  and 
overlay  roads 

-  Use  TIGER/Line 
data  and  ISL’s 
ellipsoidal  Earth 
model  toolkit 


mcarm  Data:  smi  gip  Edit  134  Detections 


20  25  30  35  40  45  50  55  60  65  70 


x  distance  (km)  ->  Longitude 


Top;  No  Edit: 
Bottom;  GIP  Edit 

Number  of 
detections  approx, 
doubled  with  GIP 
editing 

GIP  edited 
detections  appear 
more  clustered 
around  major  road 
locations 
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Summary  &  Future  Work 

•  Adaptive  thresholding  algorithm  for  multi-pass  GIP-editing  has  been 
developed  and  applied  to  both  simulated  and  experimental  data 

-  Theoretical  understanding  of  key  features  of  algorithm  developed;  robust 

-  Applicable  to  full  and  reduced  DoF  STAP  implementations 

-  Practically  implementable  with  limited  a  priori  knowledge  required 

-  Simulations  of  other  X-  and  L-band  scenarios  have  also  been  analyzed 

•  In  high-fidelity  simulations  GIP  editing  demonstrated  dramatic 
performance  improvement  over  no-editing 

-  Effectively  recover  performance  of  clutter-only  training  data 

•  Range  sidelobes  from  use  of  LFM  waveforms  degrades  GIP  editing 
performance  due  to  target  desensitization  but  significant  benefit 
over  no-editing  still  achieved  (see  associated  paper) 

•  MCARM  experimental  analysis  demonstrated  target  detection 
performance  improvement  due  to  GIP-editing 

-  Apply  to  other  experimental  data  sets,  as  available 

•  Computational  complexity  assessment  ongoing 

-  Updating  using  “sliding  window”  in  data  domain  being  examined  to  optimize 
efficiency  /  performance  trade-off 

•  Address  platform-specific  implementation  issues 
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BACK-UP  SLIDES 
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Datacube  with  LFM  Targets 

Previous  results  with  simulated  data  used  an  ideal  waveform 

Similar  analysis  performed  on  datacube  using  a  compressed 
LFM  waveform  (with  a  Taylor  nbar=6  taper)  for  all  targets  in 
the  scene  (duty  factor  =  0.1) 

Perform  sensitivity  study  to  sidelobe  level 

Compressed  Targets 


Taylor  taper,  25  dB,  nbar  =  6 


ideal  waveform 


Taylor,  25  dB,  nbar  =  6 


-50 


0  100  200  300  400  500 

range  (m) 

Produces  range 
sidelobes  from  target 
returns 

Examine  sidelobe  levels 
of  {20,  25,  30,  35,  40}  dB 


0  50  100  150 

Doppler  (m/s) 


0  50  100 

Doppler  (m/s) 
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ISL 


power  (dB) 


Pd  vs  Pfa  Performance  with  LFM 


Pfg  vs.  Pd  :  SNRq  =  25  dB,  Doppler  =  0  m/s 


Ideal  wv,  No  Edit 
Ideal  wv,  GIP  Edit 
Tayl40,  No  Edit 
Tayl40,  GIP  Edit 
Tayl35,  No  Edit 
Tayl35,  GIP  Edit 
Tayl30,  No  Edit 
Tayl30,  GIP  Edit 
Tayl25,  No  Edit 
Tayl25,  GIP  Edit 
Tayl20,  No  Edit 
Tayl20,  GIP  Edit 


Same  CFAR  analysis  as  done  previously,  targets  injected  at  Doppler  velocity 
0  m/s  (clutter  at  51  m/s) 


5-Bin  post-Doppler  with  K=110  sample  support,  0  dB  diagonal  loading 
Compare  “No  Edit”  and  “GIP  Edit”  performance  for  all  waveforms 


Both  “No  Edit”  and  “GIP  Edit”  results  degrade  with  increased  sidelobe  level 

Note  comparable  performance  for  35  and  40  dB,  consistent  with  SINR  Loss 
results 


•  Except  at  largest  sidelobe  level  (20  dB),  GIP  editing  provides  improvement  in 
performance  over  no-editing  with  LFM  waveform 
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